High resolution free-view interpolation of planar structure

ABSTRACT

A method and system is proposed to create a generative model to interpolate any view of a planar scene given a sequence of reference views and a synthesis view that is optimized by the marginalization of photometric regulation, and geometric registration parameters. According to one aspect of the claimed subject matter, a technique is proposed to combine information from varying input camera poses. Planar homography based image super resolution in free view interpolation for planar structure is applied to the combined information. Non-redundant information is combined in such a manner that the high resolution and free view problems in traditional 2D based image-based rendering techniques are overcome.

CROSS REFERENCE TO RELATED APPLICATIONS

This application claims the benefit of U.S. Provisional Application 61/916,700 to Hu et al., filed Dec. 16, 2013, which is incorporated by reference herein in its entirety and for all purposes.

BACKGROUND

In many imaging applications, users may wish to access a random viewing angle of a captured scene in high resolution, even when this specifically queried imagery is not among the set of acquired sample images. In theory, 3D reconstruction based rendering can be applied to generate such an (artificial) image. However, accurate camera calibration over large scale photo collections is needed and is highly complex in nature. Image stitching based approaches, such as panorama imaging, can also be applied. However, such schemes are unable to provide free view interpolation or resolution enhancement.

Image-based rendering (IBR) is a technique that renders novel views of the scene given a set of images, and has long been studied in computer graphics and computer vision. With the rapid progress and extensive deployment of mobile devices in recent years, there is strong consumer interest in developing light weight algorithms capable of high quality free-view interpolation with zoom-in/out effects given limited uncalibrated views.

Two widely used approaches for IBR are 3D points cloud reconstruction and image stitching, respectively. The resolution of the synthesized view of 3D reconstruction based IBR depends highly on the accuracy of the recovered points cloud. In this case, a large number of reference images from varying views are required for calibration. For image stitching based IBR, even though strict controls over camera positions are usually required, serious artifacts may still occur at the seam because of the ignorance of depth and view disparity during image registration, and due to the non-redundant information in the overlapping area of adjacent frames. Both approaches suffer from several unresolved issues when high resolution free view interpolation is desired. This is particularly true for mobile devices when limited computation resources are available and the acquisition of reference images cannot be properly controlled. In general, mobile users prefer to render a novel view of scene within a short time, based on a few available views, either taken by users themselves or acquired from the Internet.

To determine where a point appears in a specific view given its positions in a set of images, one straightforward and common way is to register all camera poses, reconstruct the 3D scene, and then project the 3D point onto the desired view. For example, where techniques for registration and rendering large scale photo collections are involved, the method may include interactively browsing and exploring large, unstructured collections of photographs. More robust 3D Six Degrees of Freedom (6DOF) registration and scene specific controls can be achieved by related techniques. The precise calibration of full 3D camera poses for all reference images is achieved by leveraging structure-from-motion algorithms with high time complexity, in which multi-core GPUs have been proposed. However, the acquisition and calibration of large scale photo collections is still a huge challenge for users with limited resource and real time needs.

Image stitching based IBR, on the other hand, usually has low time complexity, and is available even in consumer cameras, for example, to obtain a panorama shot. However, control over images is usually required for panorama generation. For standard panoramas, for example, all images are assumed to maintain the same camera center and to compose a wide field of view. On the other hand, for multi-perspective panoramas, reference images are required to correspond to a specific orbit, although the final mosaic combines views from different directions. To minimize the artifacts caused primarily by the disparities between images, dense scene depth can be estimated and structure-from-motion algorithms have also been employed for registration and determination of projection surfaces. For both types of panoramas, proper segments are cropped, registered and stitched. Unfortunately, high resolution with free view has heretofore been virtually impossible to achieve, since the redundant information in the overlapping areas of adjacent views are typically not used for anything other than seam smoothing.

When limited in planar structure, 2D based high resolution free view interpolation can be regarded as a special case of general Super Resolution (SR) problem in 3D space. Super Resolution has long been studied to generate high resolution image by combining the non-redundant information from multiple low-resolution images. Image registration and SR are often treated as distinct and sequential processes. More recently, a Bayesian approach has been developed to estimate SR by optimizing the marginal probability of the observed low-resolution images directly. Lower memory requirements and sharper results are reported by integration over the point-spread function (PSF) and motion parameters. However, all previous SR methods can only treat slight rotation and translation over reference images. These registration schemes would fail when there is disparity in either depth, pitch, or yaw in reference images.

SUMMARY

This Summary is provided to introduce a selection of concepts in a simplified form that is further described below in the Detailed Description. This Summary is not intended to identify key features or essential features of the claimed subject matter, nor is it intended to be used to limit the scope of the claimed subject matter.

As a solution to the problems noted above inherent to traditional image-based rendering techniques, the present disclosure describes a novel free view image super-resolution scheme by applying Lie group theory to 3D space to interpolate free view for planar structures. A Bayesian model is constructed for marginalizing over photometric regulation and geometric registration parameters. The proposed scheme is able to achieve desired performance against the state-of-the-art image super-resolution approaches and successfully obtain registration in full 6 degrees-of-freedom (6DOF). Compared to existing image based rendering schemes, embodiments of the present disclosure achieves free view interpolation for planar structures with higher resolution and less distortion.

According to an aspect of the present disclosure, a novel IBR scheme is described, whereby free view interpolation can be provided for planar structures given uncontrolled, limited number of input images. To obtain the position and intensity for each point in the synthesis image, instead of re-projecting from the 3D model or making copy from the original images, it is formulated as a super-resolution (SR) problem with desired registration in full 6DOF. According to one or more embodiments, a Bayesian approach is employed to determine the SR image given multiple low resolution images. To solve the challenge when significant disparities exist between reference views, the Lie group theory is applied to the Bayesian model to implement the optimization on geometric registration together with photometric regulation. With such optimization, this scheme is able to realize any zoom-in and zoom-out as well as rotation effects of the scene based on only limited views.

According to a second aspect of the present disclosure, a solution is provided for the image SR problem in full 6DOF. In contrast with general SR methods that can only treat slight rotation and translation in reference images, the scheme according to the present disclosure can handle scale change and 3D transformation effectively. According to one or more embodiments, a novel IBR scheme derived from image SR problem is disclosed. The scheme can generate free view interpolation for planar structures with higher resolution and less distortion using limited uncontrolled input images that is different from either 3D reconstruction or image stitching based IBR.

BRIEF DESCRIPTION OF THE DRAWINGS

The accompanying drawings, which are incorporated in and form a part of this specification, illustrate embodiments of the invention and, together with the description, serve to explain the principles of the invention:

FIG. 1 a depicts a diagram that illustrates the goal to generate an exemplary high resolution free view interpolation of a planar structure, in accordance with embodiments of the present invention.

FIG. 1 b illustrates multiple exemplary input views corresponding to the scene depicted in FIG. 1 a, in accordance with embodiments of the present invention.

FIG. 2 depicts a flowchart to generate an exemplary high resolution free view interpolation of a planar structure, in accordance with embodiments of the present invention.

FIG. 3 illustrates an exemplary directed graph model, in accordance with embodiments of the present invention.

FIG. 4 shows exemplary results on a text sequence, in accordance with embodiments of the present invention.

FIG. 5 depicts exemplary images of a scene panorama view interpolation on planar structures of an exemplary video sequence, in accordance with an embodiment of the present disclosure.

FIG. 6 a depicts an exemplary Jigsaw puzzle dataset, in accordance with embodiments of the present invention.

FIG. 6 b depicts exemplary results between input frames, an average image, and the image generated using high resolution free view interpolation, in accordance with embodiments of the present invention.

FIG. 7 depicts an exemplary operating environment upon which embodiments of the present disclosure may be implemented.

DESCRIPTION OF THE INVENTION

Reference will now be made in detail to several embodiments. While the subject matter will be described in conjunction with the alternative embodiments, it will be understood that they are not intended to limit the claimed subject matter to these embodiments. On the contrary, the claimed subject matter is intended to cover alternative, modifications, and equivalents, which may be included within the spirit and scope of the claimed subject matter as defined by the appended claims.

Furthermore, in the following detailed description, numerous specific details are set forth in order to provide a thorough understanding of the claimed subject matter. However, it will be recognized by one skilled in the art that embodiments may be practiced without these specific details or with equivalents thereof. In other instances, well-known processes, procedures, components, and circuits have not been described in detail as not to unnecessarily obscure aspects and features of the subject matter.

Portions of the detailed description that follow are presented and discussed in terms of a process. Although steps and sequencing thereof are disclosed in figures herein describing the operations of this process, such steps and sequencing are exemplary. Embodiments are well suited to performing various other steps or variations of the steps recited in the flowchart of the figure herein, that not all of the steps depicted may be performed, or that the steps may be performed in a sequence other than that depicted and described herein.

Some portions of the detailed description are presented in terms of procedures, steps, logic blocks, processing, and other symbolic representations of operations on data bits that can be performed on computer memory. These descriptions and representations are the means used by those skilled in the data processing arts to most effectively convey the substance of their work to others skilled in the art. A procedure, computer-executed step, logic block, process, etc., is here, and generally, conceived to be a self-consistent sequence of steps or instructions leading to a desired result. The steps are those requiring physical manipulations of physical quantities. Usually, though not necessarily, these quantities take the form of electrical or magnetic signals capable of being stored, transferred, combined, compared, and otherwise manipulated in a computer system. It has proven convenient at times, principally for reasons of common usage, to refer to these signals as bits, values, elements, symbols, characters, terms, numbers, or the like.

It should be borne in mind, however, that all of these and similar terms are to be associated with the appropriate physical quantities and are merely convenient labels applied to these quantities. Unless specifically stated otherwise as apparent from the following discussions, it is appreciated that throughout, discussions utilizing terms such as “accessing,” “writing,” “including,” “storing,” “transmitting,” “traversing,” “associating,” “identifying” or the like, refer to the action and processes of a computer system, or similar electronic computing device, that manipulates and transforms data represented as physical (electronic) quantities within the computer system's registers and memories into other data similarly represented as physical quantities within the computer system memories or registers or other such information storage, transmission or display devices.

The claimed subject matter is directed to a novel free view image super-resolution scheme by applying Lie group theory to 3D space to interpolate free view for planar structures. A Bayesian model is constructed for refining (re-calculating) photometric regulation and geometric registration parameters. The proposed scheme is able to achieve desired performance against the state-of-the-art image super-resolution approaches and successfully obtain registration in full 6 degrees-of-freedom (DOF). Compared to existing image based rendering schemes, embodiments of the present disclosure achieves free view interpolation for planar structure with higher resolution and less distortion.

According to the present disclosure, a novel IBR scheme is described, whereby free view interpolation can be provided for planar structures given uncontrolled, limited number of input images. To obtain the position and intensity for each point in the synthesis image, instead of re-projecting from the 3D model or making copy from the original images, it is formulated as a super-resolution (SR) problem with desired registration in full 6DOF. A Bayesian approach is employed to determine the SR image given multiple low resolution images. To solve the challenge when significant disparities exist between reference views, the Lie group theory is applied to the Bayesian model to implement the optimization on geometric registration together with photometric regulation parameters. With such optimization, this scheme is able to realize any zoom-in and zoom-out as well as rotation effects of the scene based on only limited views.

FIGS. 1 a, 1 b, and 2 illustrate an overview of an exemplary high resolution free view interpolation process in accordance with an embodiment of the present disclosure. FIG. 1 a depicts a diagram 100 a that illustrates the goal to generate any free view of the scene V_(f) given a set of views of a planar scene (V₁-V₁₀). FIG. 1 b illustrates multiple exemplary input views 100 b corresponding to the scene V_(f) depicted in FIG. 1 a. FIG. 2 depicts a flowchart 200 for the proposed scheme.

As depicted in FIG. 1 a, novel (or “synthesis”) view V_(f) of a planar scene can be interpolated from a given set of views (V₁-V₁₀) from different angles and depths as demonstrated in FIG. 1 b. In one or more embodiments, V_(f) can be any view such as panorama or a normal view that is not available or captured (for example, famous paintings are not always allowed to be captured from a relatively short distance in museums). The flowchart 200 of the proposed scheme is depicted in FIG. 2. As presented in FIG. 2, reference images are initialized at step 201. Initialization may include, for example, registration of the reference images by planar homography in full 6DOF and the generation of a pre-estimated image that includes a synthesis image of the reference images. In one or more embodiments, the synthesis image may be initially generated using a homography transformation of one of the reference images, and an average of color intensity values of the input images. Subsequently, the synthesis image undergoes an optimization process 203 that includes optimization of pixel coordinates and intensities of the synthesis image. Finally, new synthesis images are generated at step 205 by iteratively re-calculating geometric registration and photometric regulation parameters to refine the synthesis image. The main optimization process is based on a Bayesian generative model, which will be described below in greater detail.

In one or more embodiments, the super-resolution (SR) image is obtained by optimizing posterior distribution marginalized over rotation and translation as well as intensity parameters in the Bayesian model. A Bayesian model is adopted in the present disclosure for transformation in full 6DOF by introducing properties of the Lie group and optimizations in positions and intensity are formulated as two separate processes. With adopted registration, it is not necessary to exercise strict controls over positions of input images. Instead, it is possible to simply crop and stitch patches as traditional image stitching based IBR. Reconstructing high resolution free view is accomplished by combining information from limited views directly, in order to avoid processing large scale photo collections during 3D reconstruction.

The following sections describe in detail an embodiment of the present disclosure, which employs the proposed Bayesian model for high resolution free view interpolation. The generative model, the marginalization over photometric and regulation parameters, and the optimization procedure are presented with reference to the embodiments in detail.

Projection of the planar structure on a requested view is determined by its K available views. Assume x=(v,h) to be a set of pixels in the image of a requested view with size N, where v=[v₁, . . . , v_(N)]^(T) is an N×2 vector corresponding to coordinates of the image and h represents intensity with the dimension N×1 (which can be extended to RGB image straightforwardly). Similarly, reference pixels of x in view k is represented as y^((k))=(w^((k)),g^((k))) with w^((k))=[w₁ ^((k)), . . . , w_(N) ^((k))]^(T)(N×2) and g^((k))(N×1) representing coordinates and intensity, respectively. The relationships between x and {y^((k))} are governed by the set of planar homography matrices {H^((k))} and photometric regulation parameters {λ^((k))}.

Assuming v_(i) and w_(i) ^((k)) are coordinates of one pixel in the requested view and its reference coordinate in view k, w _(i) ^((k)) and v _(i) are used to represent the homogeneous coordinates for w_(i) ^((k)) and v_(i). With planar homography H^((k)) and Gaussian noise added to each dimension of the homogeneous coordinate, the coordinate transformation can be expressed as,

w _(i) ^((k)) =H ^((k)) v _(i)+ε_(v) ^((k))  (1)

where

$\begin{matrix} {{ɛ_{v}^{(k)} \sim {{N\left( {0,\Sigma_{1}} \right)}\Sigma_{1}}} = \begin{bmatrix} a^{- 1} & 0 & 0 \\ 0 & b^{- 1} & 0 \\ 0 & 0 & c^{- 1} \end{bmatrix}} & (2) \end{matrix}$

A linear regulation and also a Gaussian i.i.d. noise are used to represent the intensity relationships of the same scene points projected in different views. Ignoring the PSF, all reference views are of proper resolution, with the intensity g being expressed as:

g ^((k))=Δα^((k)) h+λ _(β) ^((k)) e+ε ^((k))

ε^((k)) ˜N(0,β⁻¹ I)  (3)

where e is a N×1 vector with all values equal to 1.

In one or more embodiments, the homography H^((k)) can be calculated by extracting and matching SIFT points in both images and conducting direct linear transformation (DLT) using a random sample consensus (RANSAC) method, where x′ is the pre-estimated representation of x and given x′ and y^((k)). Photometric regulation parameters {λ^((k))} can be expressed as λ^((k))=[λ_(a) ^((k))λ_(β) ^((k))]^(T), and can be calculated by linear regression given the intensities of corresponding pixels.

FIG. 3 illustrates an exemplary directed graph model 300 corresponding to the joint distribution of x′, {y^((k))}, {H^((k))}, {λ^((k))} and x. As depicted in FIG. 3, x′, the pre-estimated representation of x, is jointly distributed to photometric regular parameters λ^((k)), and homography matrix H^((k)) for each of the K available views. Reference pixels y^((k)) for a k view of the K available views is likewise jointly distributed to the photometric regulation parameters λ^((k)) and the homography matrix H^((k)), and combined with both to derive synthesis views of x, the set of pixels in the required image.

Given pre-estimated x′ and {y^((k))}, the joint distribution of geometric registration {H^((k))}, photometric regulation {λ^((k))} parameters and the new synthesis view x will be

p(x,{H ^((k)),λ^((k)) }|{y ^((k)) },x′)

p(x≦{y ^((k)),λ^((k)) ,H ^((k))})

Π_(k) p(λ^((k)) |x′,y′ ^((k)))Π_(k) p(H ^((k)) x′,y ^((k)))  (4)

To obtain the new synthesis view, {H^((k))} and {λ^((k))} are re-calculated in an integral form and the conditional distribution is calculated to obtain transformation.

$\begin{matrix} {= \begin{matrix} {p\left( {\left. x \middle| \left\{ y^{(k)} \right\} \right.,x^{\prime}} \right.} \\ {\int\mspace{14mu} {\ldots \mspace{14mu} {\int{p\left( x \middle| \left\{ {y^{(k)},\lambda^{(k)},H^{(k)}} \right\} \right)}}}} \\ {\Pi_{k}{p\left( {\left. \lambda^{(k)} \middle| x^{\prime} \right.,y^{(k)}} \right)}\Pi_{k}{p\left( {{H^{(k)}x^{\prime}},y^{(k)}} \right)}{\left\{ {H,\lambda} \right)}} \\ {\frac{p(x)}{p\left( \left\{ y^{(k)} \right\} \right.}{\int\mspace{14mu} {\ldots \mspace{14mu} {\int{\Pi_{k}{p\left( {\left. y^{(k)} \middle| x \right.,\lambda^{(k)},H^{(k)}} \right)}}}}}} \\ {{p\left( {\left. \lambda^{(k)} \middle| x^{\prime} \right.,y^{(k)}} \right)}{p\left( {\left. H^{(k)} \middle| x^{\prime} \right.,y^{(k)}} \right)}{\left\{ {H,\lambda} \right\}}} \end{matrix}} & (5) \end{matrix}$

Assuming each pixel's coordinate is independent with its intensity, Equation (5) is decomposed into multiplication of distribution for both parts.

$\begin{matrix} {p\left( {\left. x \middle| \left\{ y^{(k)} \right\} \right.,{x^{\prime}\frac{p(h)}{p\left( \left\{ g^{(k)} \right\} \right.}\Pi_{k}{\int{{p\left( {\left. g^{(k)} \middle| h \right.,\lambda^{(k)}} \right)}{p\left( {\left. \lambda^{(k)} \middle| x^{\prime} \right.,y^{(k)}} \right)}{\lambda^{(k)}}\frac{p(v)}{p\left( \left\{ w^{(k)} \right\} \right.}\Pi_{k}{\int{{p\left( {\left. w^{(k)} \middle| v \right.,H^{(k)}} \right)}{p\left( {\left. H^{(k)} \middle| x^{\prime} \right.,y^{(k)}} \right)}{H^{(k)}}}}}}}} \right.} & (6) \end{matrix}$

In one or more embodiments, the new synthesis image x is determined given x′ and {y^((k))} by searching for the intensity and the coordinates with the highest probability, respectively.

From Equation (3), for each individual view, given requested view and intensity regulation parameters, the likelihood can be expressed as:

$\begin{matrix} {{p\left( {\left. g^{(k)} \middle| h \right.,\lambda^{(k)}} \right)} \propto {\exp \left\{ {{- \frac{\beta}{2}}{\sum\limits_{n = 1}^{N}\; \left( {g_{n}^{(k)} - {\lambda_{\alpha}^{(k)}h_{n}} - \lambda_{\beta}^{(k)}} \right)^{2}}} \right\}}} & (7) \end{matrix}$

where h_(n) is the intensity value for the nth pixel.

According to one or more embodiments of the present disclosure, δ^((k)) is used to represent the difference between the real photometric regulation parameters λ^((k)) and its estimated values λ^((k)) where δ^((k))˜N(0,C) and C=y⁻¹l.

$\begin{matrix} {\begin{bmatrix} \lambda_{\alpha}^{(k)} \\ \lambda_{\beta}^{(k)} \end{bmatrix} = {\begin{bmatrix} {\overset{\_}{\lambda}}_{\alpha}^{(k)} \\ {\overset{\_}{\lambda}}_{\beta}^{(k)} \end{bmatrix} + \delta^{(k)}}} & (8) \\ {{p\left( {\left. \lambda^{(k)} \middle| x^{\prime} \right.,y^{(k)}} \right)} = {\left( \frac{C^{- 1}}{\left( {2\pi} \right)^{n}} \right)^{\frac{1}{2}}\exp \left\{ {{- \frac{1}{2}}\delta^{{(k)}T}C^{- 1}\delta^{(k)}} \right\}}} & (9) \end{matrix}$

The negative exponent term of (7) can then be a function of δ^((k)) and approximated by the second-order Taylor series.

$\begin{matrix} \begin{matrix} {{^{(k)}\left( \delta^{(k)} \right)} = {\frac{\beta}{2}{\sum\limits_{n = 1}^{N}\; \left( {g_{n}^{(k)} - {\lambda_{\alpha}^{(k)}h_{n}} - \lambda_{\beta}^{(k)}} \right)^{2}}}} \\ {= {F^{(k)} + {L^{{(k)}T}\delta^{(k)}} + {\delta^{{(k)}T}M^{(k)}\delta^{(k)}}}} \end{matrix} & (10) \end{matrix}$

In one or more embodiments, Huber prior is used for prior distribution of intensity,

$\begin{matrix} {{{p(h)} = {\frac{1}{Z_{z}}\exp \left\{ {{- \frac{v}{2}}{p\left( {{Dx},\alpha} \right)}} \right\}}}{{p\left( {z,\alpha} \right)} = \left\{ {{2\alpha {z}} - {\alpha^{2}\mspace{31mu} \begin{matrix} {{{if}{z}} < \alpha} \\ {otherwise} \end{matrix}}} \right.}} & (11) \end{matrix}$

where Dx is the gradient of x. and a and v are parameters of the Huber potential function. In one or more further embodiments, Zx is the partition function.

The photometric regulation part of (6) can then be derived as follows, assuming S^((k))=M^((k))+½C⁻¹

$\begin{matrix} {{{p(h)}{\prod\limits_{k}\; {\int{{p\left( {\left. g^{(k)} \middle| h \right.,\lambda^{(k)}} \right)}{p\left( {\left. \lambda^{(k)} \middle| x^{\prime} \right.,y^{(k)}} \right)}{\lambda^{(k)}}}}}}{\alpha \mspace{14mu} \exp \left\{ {{- \frac{v}{2}}{\rho \left( {{Dx},\alpha} \right)}} \right\}}{{\prod\limits_{k}\; {\exp \left\{ F^{(k)} \right\} {\int{\exp \left\{ {L^{{(k)}T} - {\delta^{{(k)}T}S^{(k)}\delta^{(k)}}} \right\} {\delta^{(k)}}}}}} = {{\left( {2\pi} \right)\left( \frac{N}{2} \right){S^{(k)}}} - \frac{1}{2}}}{\exp \left\{ {{- \frac{v}{2}}{\rho \left( {{Dx},\alpha} \right)}F^{{(k)}T}S^{{(k)} - 1}L^{(k)}} \right.}} & (12) \end{matrix}$

Similar to the derivations described above, the likelihood given requested view and registration parameters can be derived from (1) and (2):

$\begin{matrix} {{p\left( {\left. w^{(k)} \middle| v \right.,H^{(k)}} \right)}\alpha \; \exp \left\{ {1 - {\frac{1}{2}{\sum\limits_{n = 1}^{N}\; {\left( {{\underset{\_}{w}}_{n}^{(k)} - {H^{(k)}{\underset{\_}{v}}_{n}}} \right)^{T}{\sum\limits_{1}\; {{- 1}\left( {{\underset{\_}{w}}_{n}^{(k)} - {H^{(k)}{\underset{\_}{v}}_{n}}} \right)}}}}}} \right\}} & (13) \end{matrix}$

Since the planar homography H(k)εSL(3)—which is a space of 3×3 real matrices with 8-DOF belonging to the Lie group—the difference between the real registration matrix H(k) and its estimation H(k) can be expressed as follows:

$\begin{matrix} {{{{{\overset{\_}{H}}^{{(k)} - 1}H^{(k)}} = {\text{?} = {{1 + {\text{?}\frac{1}{n!}\text{?}}} \approx {1 + G^{(k)}}}}}{\text{?}\text{indicates text missing or illegible when filed}}}} & (14) \end{matrix}$

with the Gaussian distribution of H^((k)) on SL(3) is defined as:

$\begin{matrix} {{{{p\left( H^{(k)} \middle| {\text{?}y^{(k)}} \right)} \times \exp \left\{ {{- \frac{1}{2}}f\text{?}\Sigma_{2}^{- 1}{f\left( G^{(k)} \right)}} \right\}}{\text{?}\text{indicates text missing or illegible when filed}}}} & (15) \end{matrix}$

where f: sl(3)→

⁸ is a mapping that projects the tangent space of SL(3) to an 8 dimensional space and Σ₂=d⁻¹l is a 8×8 diagonal matrix.

In one or more embodiments, assuming X^((k))=f(G^((k))) is an 8×1 vector, the negative exponent term of (13) can be approximated by a second-order Taylor series over X^((k)) expressed as:

$\begin{matrix} {{\begin{matrix} {{\text{?}\left( X^{(k)} \right)} = {\frac{1}{2}\text{?}\left( {\text{?} - \text{?}} \right){\Sigma_{1}^{- 1}\left( {\text{?} - {H^{(k)}\text{?}}} \right)}}} \\ {= {F^{(k)} + {\text{?}{\overset{\_}{M}}^{(k)}X^{(k)}}}} \end{matrix}{\text{?}\text{indicates text missing or illegible when filed}}}\mspace{284mu}} & (16) \end{matrix}$

Since the coordinates of the synthesis view are uniformly distributed, p(v) can be considered as constant.

The geometric registration part of (6) can then be derived as follows, assuming Ŝ^((k))=ij^((k))+⅓Σ₂ ⁻¹:

$\begin{matrix} {{{{P(v)}{\prod\limits_{k}\; {\int{{p\left( w^{(k)} \middle| {v \cdot H^{(k)}} \right)}{p\left( H^{(k)} \middle| {x^{\prime} \cdot y^{(k)}} \right)}{H^{(k)}} \times \left( {2\pi} \right)\text{?}{{\overset{\_}{S}}^{(k)}}^{\frac{1}{2}}\exp \left\{ {{\overset{\_}{F}}^{(k)}\_ \text{?}{\overset{\_}{S}}^{{(k)} - 1}{\overset{\_}{L}}^{(k)}} \right\}}}}}{\text{?}\text{indicates text missing or illegible when filed}}}\mspace{301mu}} & (17) \end{matrix}$

letting

be the objective function as the log of (6) and neglecting the constant terms leaves:

$\begin{matrix} {{{\mathcal{L} = {{\text{?}{\rho \left( {{Dx} \cdot n} \right)}} + {\text{?}\left( {F^{(k)} + {\frac{1}{2}\log {S^{(k)}}} - {\text{?}S^{{(k)} - 1}L^{(k)}} + {\overset{\_}{F}}^{(k)} + {\frac{1}{2}\log {{\overset{\_}{S}}^{(k)}}} - {\text{?}S^{{(k)} - 1}{\overset{\_}{L}}^{(k)}}} \right)}}}{\text{?}\text{indicates text missing or illegible when filed}}}\mspace{290mu}} & (18) \end{matrix}$

This objective function can be optimized using bounded quasi-Newton method of BFGS (L-BFGS-B). The gradients of

over h and v can be expressed as

$\begin{matrix} {\text{?} = {{\text{?}\frac{}{{Dx}}{\rho ({Dx})}} + {\sum\limits_{k = 1}^{K}\; \left\{ {\text{?} - {\frac{1}{2}\text{?}S^{{(k)} - 1}\text{?}} + {\left( {{\frac{1}{2}\text{?}} + {\text{?}\left( {\text{?}S^{{(k)} - 1}\text{?}S^{{(k)} - 1}} \right)}} \right)\text{?}}} \right\}}}} & (19) \\ {{{\text{?}{\sum\limits_{k = 1}^{K}\; \left\{ {\text{?} - {\frac{1}{2}\text{?}{\overset{\_}{S}}^{{(k)} - 1}\text{?}} + {\left( {{\frac{1}{2}\text{?}} + \left( {\text{?}{\overset{\_}{S}}^{{(k)} - 1}\text{?}{\overset{\_}{S}}^{{(k)} - 1}} \right)} \right)\text{?}}} \right\}}}{\text{?}\text{indicates text missing or illegible when filed}}}\mspace{301mu}} & (20) \end{matrix}$

F, L, S, F, L, S and their gradients can be obtained analytically.

The process above is depicted in the flowchart of FIG. 2. Initially, the planar homography between all reference views is calculated and a synthesis view is generated at step 201. In one or more embodiments, the synthesis view may be generated based on photometric regulation and geometric registration parameters derived from reference views (input images). Geometric registration parameters (including, for example, one or more planar homography transformation matrices between the reference views) may be calculated by using SIFT matching and assuming the synthesis view is transformed from one reference view. Thereafter, the planar homography between the synthesis view and all reference views is obtained at step 203 to optimize the pixel intensities and coordinates of the synthesis view. In one embodiment, photometric regulation parameters include color and intensity values for pixels in the synthesis image. The initial intensity for each pixel in the synthesis view may be determined based on the average of all reference pixels. Subsequently, the initial photometric regulation parameters and geometric registration parameters are iteratively recalculated at step 205 to refine the generated synthesis image. In one or more embodiments, the photometric regulation parameters may be re-calculated using linear regression.

When the synthesis image is large and each pixel can only find reference pixels in partial of reference views, the performance may be influenced if all pixels of the image are considered at one time. In the presented embodiments, an image for the requested view is synthesized slit by slit. That is to say, positions and intensity for pixels in the same column may be optimized at the same time. After all pixels of the image are refined geometric registration {H^((k))} and photometric regulation {λ^((k))} parameters are recalculated, and a new synthesis image is generated based on the re-calculated parameters

Exemplary results obtained from application of an exemplary high resolution free view interpolation process in accordance with embodiments of the present disclosure on three different sample datasets are presented herein: text for image SR reconstruction, a self-collected jigsaw puzzle dataset and a planar structure of a video sequence for long scene panorama. Intensity and coordinate values are normalized in the experiments. In one or more embodiments, the results obtained are insensitive to the parameter setting. The parameters that perform well on all datasets are set as follows: the Huber parameter a is fixed at 0.01 and v is set to 1. Intensity values may differ a lot from image to image because of the distinction in lighting conditions, so a larger variance is selected for their distribution comparing with registration parameters. The variance for the distribution of intensity value β and γ are set to 5 and the variance for the distribution of registration parameters a, b, c and d are set to 0.1. As demonstrated for experiments on image super-resolution dataset, the presented scheme not only can compete with the state-of-the-art SR methods but also can achieve image registration in full 6DOF.

FIGS. 4( a)-(g) shows exemplary results 400 on a text sequence from MDSP SR and demosaicing datasets, which contains 30 uncompressed low resolution grayscale frames of size 57×49. As depicted in FIG. 4, sequence 400 includes (a) 1/30 low resolution frame; (b) average of input images after feature matching based registration; (c) bi-cubic interpolation image; (d) result by Bayesian SR; (e) results by the proposed scheme; (f) result by Bayesian SR when enrolling scale change and 3D rotation; (g) results by the described scheme when enrolling scale change and 3D rotation.

In one or more embodiments, FIG. 4( a) may be one of these 30 input frames. FIGS. 4( b)-(g) present results by increasing the resolution with a factor 2. To set up the baseline, the homography is calculated between frames through SIFT feature matching and the average all input frames after registration as FIG. 4( b). The bicubic interpolation image (depicted as FIG. 4( c) is calculated as well.

The comparison between the Bayesian super-resolution and an embodiment of the present disclosure on the text sequence are depicted in FIG. 4( d) and FIG. 4( e). The proposed scheme exhibits competent sharp results that compares favorably to the Bayesian SR scheme. Another group of experiments are carried out by enrolling scale change and rotation along all 3 axes in the original text sequence. The results from a comparison between Bayesian SR in accordance with the prior art and the present disclosure is shown as FIG. 4( f) and FIG. 4( g). Marginalizing over 3D planar homography, the scheme according to the present disclosure can still obtain sharp results, while the results by Bayesian SR become blurred since only 2D rotation and translation is considered.

FIG. 5 illustrates an resulting image 500 of scene panorama view interpolation on planar structures of an exemplary video sequence in accordance with an embodiment of the present disclosure. As presented in FIG. 5, the video sequence consists of 1800 frames with each frame being 360×240 pixels in size. Selected frames from the video sequence are extracted and the results depicted in FIG. 5. The results of image-based rendering with high resolution free view interpolation methods as described herein are depicted in FIGS. 5( d)(e)(f), juxtaposed with results from minimal aspect distortion (MAD) mosaicing methods (depicted in FIGS. 5( a)(b)(c)). MAD mosaicing uses dense depth to minimize the geometrical distortions, which gives excellent performance for long scene with varying depths. However, it is clear that the distortion at the seam between two cropped patches is still hard to avoid, even though the two cropped patches belong to the same plane. The proposed scheme described herein successfully overcomes the distortion at the frame of the door, comparing results in FIG. 5( b) and FIG. 5( e). FIG. 5( c) and FIG. 5( f) depict an inside view of the room. FIG. 5( g) and FIG. 5( h) present views with focus axis rotated 10° along the vertical and horizontal axes from perpendicular to the wall. Since the long panorama by MAD is the pure combination of patches of input images, higher resolution as well as view change cannot be achieved.

FIG. 6 a illustrates an exemplary Jigsaw puzzle dataset 601 a containing twenty four (24) 400×300 images of a completed jigsaw puzzle which is hung on a white wall in accordance with an embodiment of the present disclosure. The 24 images are captured from varying angles and depths, that may contain incomplete information required for novel view synthesizing. As depicted in FIG. 6 a, the positions and images of certain viewpoints are similar to the positions and images of corresponding viewpoints described above with respect to FIG. 1. The interpolated view V_(f) with size 800×600 by the proposed scheme is shown in FIG. 6 a with the requested pose, correct intensity and clear edge. The comparison results on enlarged details are shown in FIG. 6 b. Selected patches of input views are presented in the first row 601 b of FIG. 6 b, exhibiting low resolution and noise such as reflection or shadow. For a baseline, the average image is obtained by registering over all input images using SIFT matching. Blurring effects and double edges of jigsaw puzzle pieces are shown in the second row 603 b of FIG. 6 b due to registration errors. The results are shown in the third row 605 b of FIG. 6 b, which shows clear piece edges and correct intensity without disturbance by reflection, benefiting from marginalization over planar homography and photometric regulation.

Embodiments of the claimed subject matter are directed to methods and systems for image rendering using a high resolution free view interpolation method. FIG. 7 shows an exemplary operating environment in accordance with one embodiment of the present invention. Digital system 700 depicts the components of a basic system in accordance with embodiments of the present invention providing the execution platform for certain hardware-based and software-based functionality. Although specific components are disclosed in system 700 it should be appreciated that such components are examples. That is, embodiments of the present invention are well suited to having various other components or variations of the components recited in system 700. It is appreciated that the components in system 700 may operate with other components other than those presented, and that not all of the components of system 700 may be required to achieve the goals of system 700.

In a typical embodiment, System 700 includes sensor 703, image signal processor (ISP) 705, memory 707, input module 709, central processing unit (CPU) 711, display 713, communications bus 715, and power source 716. Power source 716 supplies power to system 700 and may, for example, be a DC or AC power source. CPU 711 and the ISP 705 can also be integrated into a single integrated circuit die and CPU 711 and ISP 705 may share various resources, such as instruction logic, buffers, functional units and so on, or separate resources may be provided for image processing and general-purpose operations. System 700 can be implemented as, for example, a digital camera, cell phone camera, portable device (e.g., audio device, entertainment device, handheld device), webcam, video device (e.g., camcorder) and the like.

Sensor 703 receives light via a lens 701 and converts the light received into a signal (e.g., digital or analog). According to some embodiments, lens 701 may be permanently attached to the system 700. Alternatively, lens 701 may be detachable and interchangeable with lens of other properties. These properties may include, for example, focal lengths, apertures and classifications. In typical embodiments, lens 701 may be constructed of glass, though alternate materials such as quartz or molded plastics may also be used. Sensor 703 may be any of a variety of optical sensors including, but not limited to, complementary metal-oxide-semiconductor (CMOS) or charge-coupled device (CCD) sensors. Sensor 703 is coupled to communications bus 715 and may provide image data received over communications bus 715. In further embodiments, sensor 703 includes light intensity sensing capability, and the image data received may include data corresponding to the determined intensity of the light in a scene or image.

Image signal processor (ISP) 705 is coupled to communications bus 715 and processes the data generated by sensor 703. More specifically, image signal processor 705 processes data from sensor 702 for storing in memory 707. For example, image signal processor 705 may compress and determine a file format for an image to be stored in within memory 707.

The input module 709 allows the entry of user-input into system 700 which may then, among other things, control the sampling of data by sensor 703 and subsequent processing by ISP 705. Input module 709 may include, but it not limited to, navigation pads, keyboards (e.g., QWERTY), buttons, touch screen controls (e.g., via display 113) and the like.

The central processing unit (CPU) 711 receives commands via input module 709 and may control a variety of operations including, but not limited to, sampling and configuration of sensor 703, processing by ISP 705, and management (e.g., the addition, transfer, and removal) of images and/or video from memory 707. In one embodiment, reference images 717 and programmed instructions may be stored in memory 707. CPU 711 may be configured to execute the plurality of programmed instructions to render a high resolution free view image of said captured scene. For example, CPU 711 may be configured to: receive a requested view point; generate an initial synthesis image based on the plurality of input images 717 in the memory 707; calculate a planar homography corresponding to said plurality of reference view points; optimize said initial synthesis image as a new synthesis image representing a planar image of said captured scene in said request view point based on an interpolation of said plurality of input images and said planar homography; and store the new synthesis image back in memory 707.

According to the present disclosure, a solution has been described for efficient generation of a novel image in high resolution and in full 6DOF. In contrast with general Super Resolution methods that can only treat slight rotation and translation in reference images, the scheme according to the present disclosure can handle scale change and 3D transformation effectively. Furthermore, a novel IBR scheme derived from image SR problem has been disclosed. The scheme can generate free view interpolation for planar structures with higher resolution and less distortion using limited uncontrolled input images that is different from either 3D reconstruction or image stitching based IBR.

Although certain preferred embodiments and methods have been disclosed herein, it will be apparent from the foregoing disclosure to those skilled in the art that variations and modifications of such embodiments and methods may be made without departing from the spirit and scope of the invention. It is intended that the invention shall be limited only to the extent required by the appended claims and the rules and principles of applicable law. 

What is claimed is:
 1. A computer implemented method of rendering a free view image of a captured scene, said method comprising: accessing a plurality of input images stored in a memory of a computing device and corresponding to a plurality of reference view points, wherein an input image represents a planar image of said captured scene at a respective reference view point; receiving a requested view point; generating, by a processor of the computing device, an initial synthesis image based on said plurality of input images; calculating, by the processor of the computing device, geometric registration parameters and photometric regulation parameters corresponding to said plurality of input images; and generating, by the processor of the computing device, a new synthesis image representing a planar image of said captured scene in said request view point based on an interpolation of said plurality of input images, said geometric registration parameters, and said photometric regulation parameters.
 2. The computer implemented method of claim 1, wherein calculating said geometric registration parameters comprises determining a positional relativity between pixels comprised in said initial synthesis image and pixels comprised in said plurality of input images.
 3. The computer implemented method of claim 2, wherein said determining said positional relativity comprises generating a planar homography transformation matrix relating the plurality of input images.
 4. The computer implemented method of claim 1, wherein said photometric regulation parameters corresponds to a relationship between intensity values of pixels comprised in said initial synthesized image with intensity values of pixels comprised in said plurality of input images.
 5. The computer implemented method of claim 1, wherein said new synthesis image comprises initial photometric regulation parameters, said initial photometric regulation parameters being calculated based on average intensity values of corresponding pixels of said plurality of input images.
 6. The computer implemented method of claim 1, further comprising: iteratively re-calculating said photometric regulation parameters based on said new synthesis image; iteratively re-calculating said geometric registration parameters based on said new synthesis image; and re-estimating said new synthesis image based on said re-calculated photometric regulation parameters and said re-calculated geometric registration parameters.
 7. The computer implemented method of claim 6, wherein said iteratively re-calculating said photometric regulation parameters comprises iteratively re-calculating said photometric regulation parameters in accordance with a Bayesian method.
 8. The computer implemented method of claim 6, wherein said iteratively re-calculating said geometric registration parameters comprises iteratively re-calculating said geometric registration parameters in accordance with a Lie group method.
 9. The computer implemented method of claim 5, wherein said re-calculating said photometric regulation parameters and said re-calculating said geometric registration parameters of pixels are performed substantially contemporaneously
 10. The computer implemented method of claim 1, wherein said free view image has a greater resolution than each input image; and wherein said requested view point is substantially different from said reference view points associated with said plurality of input images in full 6 degrees of freedom (6DOF).
 11. The computer implemented method of claim 10, wherein said geometric registration parameters comprises a full 6DOF planar homography transformation matrix.
 12. The computer implemented method of claim 1, wherein said generating said set of geometric registration parameters and said photometric regulation parameters comprises: calculating planar homographic transformation matrices by extracting and matching SIFT points in said initial synthesis image and said plurality of input images in accordance with direct linear transformation (SLT) using a RANSAC method; and calculating photometric regulation parameters in accordance with a linear regression method.
 13. The computer implemented method of claim 12, wherein said planar homographic transformation matrices comprise 3×3 real matrices with unit determinant and 8-DOF corresponding to a Lie group.
 14. The computer implemented method of claim 12 further comprising: calculating a joint probability distribution of said planar homographic transformation matrices and said photometric regulation parameters and said initial synthesis image, wherein said generating said initial synthesis image further comprises searching said photometric regulation parameters and said geometric registration parameters based on said joint probability distribution.
 15. The computer implemented method of claim 14 further comprising: calculating an object function of said joint probability distribution; and optimizing said object function in accordance with a bounded quasi-Newton method.
 16. The computer implemented method of claim 1 wherein calculating said geometric registration parameters comprises: calculating an initial planar homography transformation matrix between said plurality of input images using SIFT matching; generating said initial synthesis image using said homography transformation matrix over a selected input image; and deriving an new planar homography transformation matrix between said initial synthesis image and said plurality of input images.
 17. The computer implemented method of claim 1, wherein said reference view points associated with said plurality of input images are substantially different in view angles and view depths.
 18. The computer implemented method of claim 1, wherein said calculating said geometric registration parameters comprises adding Gaussian i.i.d. noise to each dimension of a homogeneous coordinate.
 19. The computer implemented method of claim 1, wherein said plurality of input images comprises at least one image of the group of images consisting of: an image captured by an imaging device; a computer simulated image; a cropped portion of an image; a segmented portion of an image; an image extracted from a video file.
 20. The computer implemented method of claim 1, wherein said generating said synthesis image comprises generating said synthesis image in accordance with a super resolution (SR) method.
 21. The computer implemented method of claim 1, wherein said free view image corresponds to at least one of: a zoomed-in image, a zoomed-out image, and a rotated image of said captured scene.
 22. A system for rendering a high resolution free view image of a captured scene in an image capture device, said system comprising: an image capture device comprising a lens and a plurality of shutters, said image capture device being configured to generate a plurality of input images corresponding to a plurality of reference view points, wherein an input image represents a planar image of said captured scene at a respective reference view point; a memory, communicatively coupled to said image capture device, and configured to store image data corresponding to said plurality of input images and a plurality of programmed instructions; and a processor, communicatively coupled to the memory, and configured to execute the plurality of programmed instructions to render a high resolution free view image of said captured scene, wherein the processor is configured to receive a requested view point; generate an initial synthesis image based on the plurality of input images in the memory; calculate geometric registration parameters and photometric regulation parameters corresponding to said plurality of input images; and generate a new synthesis image representing a planar image of said captured scene in said request view point based on an interpolation of said initial synthesis image, said geometric registration parameters, and said photometric regulation parameters, wherein said geometric registration parameters correspond to a positional relativity between pixels comprised in said initial synthesis image and pixels comprised in said plurality of input images, further wherein said photometric regulation parameters corresponds to a relationship between intensity values of pixels comprised in said intensity synthesized image with intensity values of pixels comprised in said plurality of input images.
 23. The image capture device of claim 22, wherein said processor is further configured to: re-calculate said photometric regulation parameters of said plurality of reference view points based on said new synthesis image; and re-calculate said geometric registration parameters of said plurality of reference view points based on said new synthesis image.
 24. The image capture device of claim 23, wherein said geometric registration parameters comprise a planar homography transformation matrix relating the plurality of input images
 25. The image capture device of claim 23, wherein the processor is further configured to optimize said photometric regulation parameters, and said geometric registration parameters to generate said free view image representing a planar image of said captured scene in said requested view point.
 26. A non-transitory computer readable medium containing program instructions embodied therein for rendering a high resolution free view image of a captured scene, said program instructions comprising: instructions to access a plurality of input images corresponding to a plurality of reference view points, wherein an input image represents a planar image of said captured scene at a respective reference view point; instructions to receive a requested view point; instructions to generate an initial synthesis image based on said plurality of input images; instructions to calculate geometric registration parameters and photometric regulation parameters corresponding to said plurality of reference view points; and instructions to generate a new synthesis image representing a planar image of said captured scene in said request view point based on an interpolation of said plurality of input images, said geometric registration parameters, and said photometric regulation parameters.
 27. The non-transitory computer readable medium of claim 26, wherein said instructions to calculate geometric registration parameters and photometric regulation parameters comprise: instructions to calculate a relationship between intensity values of pixels comprised in said intensity synthesized image with intensity values of pixels comprised in said plurality of input images; and instructions to calculate a positional relativity between pixels comprised in said initial synthesis image and pixels comprised in said plurality of input images.
 28. The non-transitory computer readable medium of claim 27, further comprising instructions to optimize said new synthesis image, said instructions comprising: instructions to iteratively re-calculate said photometric regulation parameters corresponding to said new synthesis image; and instructions to iteratively re-calculate said geometric registration parameters corresponding to said new synthesis image.
 29. The non-transitory computer readable medium of claim 28, further comprising instructions to repeat said instructions to optimize said new synthesis image to generate said free view image representing a planar image of said captured scene in said requested view point. 